Renormalization Group Approach for the Wave Packet Dynamics in Golden-Mean 

and Silver-Mean Labyrinth Tilings 
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We study the quantum diffusion in quasiperiodic tight-binding models in one, two, and three 
dimensions. First, we investigate a class of one-dimensional quasiperiodic chains, in which the 
atoms are coupled by weak and strong bonds aligned according to the metallic-mean sequences. 
The associated generalized labyrinth tilings in d dimensions are then constructed from the direct 
product of d such chains, which allows us to consider rather large systems numerically. The electronic 
transport is studied by computing the scaling behavior of the mean square displacement of the wave 
packets with respect to time. The results reveal the occurrence of anomalous diffusion in these 
systems. By extending a renormalization group approach, originally proposed for the golden-mean 
chain, we show also for the silver-mean chain as well as for the higher-dimensional labyrinth tilings 
that in the regime of strong quasiperiodic modulation the wave-packet dynamics are governed by 
the underlying quasiperiodic structure. 

PACS numbers: 71. 23. Ft, 71.15.-m, 72.15.-v 



I. INTRODUCTION 

Understanding the relations between the atomic struc- 
ture and the physical properties of materials remains one 
of the elementary questions of condensed- matter physics. 
One research line in this quest started with the discov- 
ery of quasicrystals by Shechtman et al. in 1982»i It 
soon became clear that quasicrystals with their 5-, 8-, 
10- or 12-fold rotational symmetries, forbidden according 
to classical crystallography, can be described in terms of 
mathematical models of aperiodic tilings of a plane pro- 
posed by Penrose and Ammann in the 1970si^ Further, 
quasicrystals are often regarded to have a degree of order 
intermediate between crystals and disordered systems, 
which can be understood by considering the repetitive- 
ness of local patterns. While any local pattern in crystals 
repeats for integer combinations of the lattice vectors, 
Conway's theorem states that in a quasicrystal a given 
local pattern in a region of some diameter L will be re- 
peated within a distance of two diameters 2LA In liquids 
or amorphous systems the distance of identical local pat- 
terns of size L grows exponentially with L. Moreover, 
the deterministic construction rules of the quasicrystal 
models can also be interpreted as some kind of order in 
contrast to the situation in disordered systems. On the 
other hand, translational symmetry is missing in contrast 
to crystals. 

Various experimental investigations revealed rather ex- 
otic electrical, magnetic, and optical properties of these 
materialsi^i^ As many quasicrystals contain a high per- 
centage of well-conducting elements, it was originally as- 
sumed that their thermal and electrical transport prop- 
erties were similar to those in crystalline or amorphous 
metals However, it turned out that they rather behave 
like semiconducting or insulating materials because they 
often possess a large resistivity for small temperatures 
and the resistivity decreases with increasing temperature 
and increasing structural order of the materialsi^"— 



This motivated extensive research to get a better theo- 
retical understanding of the structure and physical prop- 
erties of quasicrystals. Today several exact results are 
known for one-dimensional quasiperiodic systemsJ^"— 
However, the properties in two or three dimensions have 
been clarified to much lesser degree and are mainly based 
on numerical studies for limited systems with a few thou- 
sands sites due to the missing translational symmetry of 
quasicrystals. For instance, the electronic properties were 
determined numerically for the Penrose tiling, the oc- 
tagonal tiling, and the Ammann-Kramer-Neri tilingJ^ii^ 
However, this makes it difficult to understand the nature 
of eigenstates and properties of macroscopic systems. To 
address this challenge, different approaches have been fol- 
lowed. For example, by expanding the time-evolution op- 
erator in Chebyshev polynomials the spreading of energy- 
filtered wave packets was studied for generalized Rauzy 
tilings with up to 10^ sitesJ^ii^ However, this approach 
only works for systems with a rather smooth local density 
of states. Further, it is possible to calculate some exact 
eigenstates of the tight-binding Hamiltonian on the Pen- 
rose tiling.— 

We follow another approach by studying the wave- 
packet dynamics for d-dimensional quasicrystalline mod- 
els with a separable Hamiltonian in a tight-binding 
approachji^ This method is based on mathematical se- 
quences, constructed by an inflation rule V describing 
the weak and strong couplings of atoms in a quasiperi- 
odic chain. The higher-dimensional labyrinth tilings are 
then constructed as a direct product of these chains and 
their eigenstates are directly calculated by multiplying 
the energies or wave functions of these chains. This al- 
lows us to study very large systems in higher dimensions 
with up to 10^° sites based on the solutions in one di- 
mension. However, the property of separability comes at 
the price that the labyrinth tilings do not possess a non- 
crystallographic rotational symmetry. Nevertheless, the 
labyrinth tiling is able to describe a subset of points of 
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the octagonal quasicrystal associated to the structure of 
actual quasicrystals with 8-fold rotational symmetry.— 

In this paper we relate the hierarchical properties of 
these quasiperiodic systems to their electronic transport 
properties by numerical calculations and a rcnormal- 
ization group (RG) approach in the regime of strong 
quasiperiodic modulation of the bond strength. Up 
to now this method has only been used to describe 
the quantum diffusion for the one-dimensional Fibonacci 
chain,—"— whereas in this paper we will also study the 
wave-packet dynamics for the silver-mean chain and the 
generalized labyrinth tilings. 

The upcoming sections are structured as follows: 
In Sec. [IT] we introduce the metallic-mean chains and 
describe the construction of the higher-dimensional 
labyrinth tilings. In Sec. IIIII wc then present numerical 
results for wave-packet dynamics in these quasiperiodic 
tilings. Using an RG approach and perturbation theory 
we derive also analytical results for the scaling behavior 
of the width of the wave packet in Sees. IIVI to IVIi which 
show a good agreement with the numerical results in the 
regime of strong quasiperiodic modulation. The results 
are briefly summarized in Sec. IVIII 
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FIG. 1. Golden-mean labyrinth tiling C^'^ constructed from 
two golden- mean chains C^", which are perpendicular to 
each other. The bond strengths of the labyrinth tiling equal 
the product of the associated bond strengths of the one- 
dimensional chains as denoted for the three different bond 
types in the left lower corner. 



II. GENERALIZED LABYRINTH TILINGS 

The construction of the d-dimensional generalized 
labyrinth tilings is based on d so-called metallic-mean 
quasiperiodic sequences, which for a parameter b are de- 
fined by the inflation rule^i 
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Starting with the symbol s we obtain after a iterations 
the ath order approximant Ca of the quasiperiodic chain. 
The length fa of an approximant Ca is given by the recur- 
sive rule fa = bfa-i+fa-2 with /o = /i = 1. Further, the 
ratio of the lengths of two successive iterants as well as 
the ratio of the numbers # of occurrence of the symbols 
w and s in an approximant approach different metallic 
means for a oo^ Thus, depending on the parameter 
b we obtain with the continued fraction representation 
Tb — \b] = [b,b,b, ...] the relations 
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fa ^^"^^ 

fa-1 #s (Ca) 



n 
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In this paper we only consider the cases 6=1 and 6 = 2. 
Thereby, the first case corresponds to the well known 
Fibonacci sequence with the golden mean tau = [1] = 
(1 + V5)/2, which is related to real quasicrystals with 5- 
and 10-fold symmetry;^ and the latter case results in the 
octonacci sequence with silver mean TAg = [2] = 1 + ^/2, 
which is related to octagonal quasicrystals <^ 

Our model describes an electron hopping from one ver- 
tex of a quasiperiodic chain to a neighboring one. The 



aperiodicity is given by the underlying quasiperiodic se- 
quence of couplings, where the symbols w and s denote 
the weak and strong bonds in the chain. Solving the 
corresponding time-independent Schrodingcr equation 

H I*') = E' I*') ^ ti_,j'fl_, + = E'^ (3) 

for the quasiperiodic systems with zero on-site potentials, 
we obtain discrete energy values E'' and wave functions 
I 'I'*) = 10 represented in the orthogonal basis 

states |/) associated to a vertex I. The hopping strength t 
in the Schrodingcr equation is given by the quasiperiodic 
sequence Ca with ts = s for a strong bond and tw = w for 
a weak bond (0 < w < s)J^^ Applying free boundary 
conditions the number of vertices is Na = /a + 1- 

The results show that the eigenvalues are symmetric 
with respect to 0. Thus, for even system sizes Na all 
energies E have a symmetric counterpart —E, but for 
odd Na there is one state E = Q, which has no corre- 
sponding state. Additionally, the eigenfunctions possess 
a symmetryi^l The eigenstate ^' with the eigenvalue E 
and the eigenstate 4* with the respective eigenvalue —E 
only differ by an alternating sign depending on the ver- 
tex I according to ^I*; = (— 1)'^*;. Further, for odd Na the 
eigenvector 4'; associated to the eigenvalue E = has a 
special structure, i.e., vanishes on all even sites 

In a next step the generalized labyrinth tilings in 
d dimensions are constructed from the Euclidian product 
Ca X Ca X ... X Ca of d quasipcriodic chains Ca r^t^ From 
a geometrical point of view these chains are perpendicu- 
lar to each other and the diagonals of the resulting grid 
correspond to the bonds of the labyrinth tiling as shown 
in Fig. [TJ Due to this product approach also the bond 
strengths of the labyrinth tiling equal the products of 
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the corresponding bond strengths of the one-dimensional 
chains as visualized in Fig. [1] Hence, in d dimension d-{-l 
different couplings occur. Note that the grid decomposes 
into 2'^"^ separate grids depending on the starting point. 
Each of these grids corresponds to a finite ath order ap- 
proximant Ca of the generalized labyrinth tiling 

The Hamiltonian of the higher-dimensional generalized 
labyrinth tiling is separable. Thus, the eigenstates of the 
labyrinth in d dimensions can be constructed from the 
product of the eigenstates of d one-dimensional chainsi^l 



puting the mean square displacement (also called width) 



dk 



^lm...n '-^ ■ ■ ^ 



dl 



(4a) 
(4b) 



The superscripts s ~ {i,j^...,k) enumerate the eigen- 
values E and r = (l^m, . . . ,n) represent the coordi- 
nates of the vertices in the quasiperiodic tiling. Due to 
the symmetries of the eigenstates of the one-dimensional 
chains mentioned above, some of the wave functions $^ 
and related eigenvalues E^ are identical and we have to 
choose a subset which only includes distinct eigenstates. 
The construction of such a subset is described in detail 
elsewhere i^i^ 



III. MEAN-SQUARE DISPLACEMENT OF A 
WAVE PACKET 



As outlined in the introduction, the connection of the 
transport properties and the quasiperiodic structure of a 
system is not yet fully understood, and thus the investiga- 
tion of transport properties of these materials continues 
to be of special interest. In this section, we study the 
wave-packet dynamics in the metallic-mean systems, in 
particular in the limit of weak coupling (w <C s). The 
results are then related to the structure of the systems 
by extending an RG approach in the following section. 

We investigate the transport properties of the 
quasiperiodic systems by means of the time evolution of a 
wave packet |T(ro, t)) = J2rec T~r(ro7 which is con- 
structed from the solutions $r(i) of the time-dependent 
Schrodinger equation H(r,t)$r(i) — ift$r(i)- The time- 
dependent wave functions can be easily obtained via 
the separation approach from the solutions of the time- 
independent Schrodinger equation in Eq. (|4|) according 
to ^^{t) = $^e^'^ *. The wave packet is then repre- 
sented as a superposition of these orthonormal eigen- 
states where we assume that it is initially lo- 
calized at the position Tq of the quasiperiodic tiling, i.e.. 



Tr(ro,t = 0) 



Hence, with the completeness rela- 



tion and normalized basis states the wave packet is de- 
fined by 



Tr(ro,t) - 



(5) 



Besides calculating the expansion of the wave packet in 
space, a more detailed analysis can be obtained by com- 



d(ro,<) 



Y,\r-ro\'\TAt)\' 



(6) 



of the wave packet. The spreading of the width d{rQ,t) 
of the wave packet in an infinite system shows anomalous 
diffusion for t — >• oo according to d{rQ,t) oc t^^'^°\ where 
the scaling exponent /3(ro) depends on the initial posi- 
tion ro of the wave packet.— '^'^'^'^ The properties of 
the electronic transport are governed by the wave-packet 
dynamics averaged over different initial positions, that is, 
we determine 



d{t) = (d(ro,0) (X 



(7) 



Thereby, the scaling exponent /? is related to the con- 
ductivity a via the generalized Drude formula for zero- 
frequency conductivity^ ' 33,34 



(8) 



The quantity e denotes the elementary charge of an elec- 
tron, q{E-p) the density of states at the Fermi level E-p. c 
a constant, and igc a characteristic time beyond which 
propagation becomes diffusive due to scatteringi^i^^ 
Hence, /3 = corresponds to the absence of diffusion, /3 = 
1/2 to classical diffusion, and /3 = 1 to ballistic spreading. 
For quasiperiodic structures one often observes anoma- 
lous diffusion characterized by < /3 < 1.— liiiil 

Within the range of anomalous diffusion, the sub- 
diffusive regime for /? < ^ is most interesting because 
in this case the conductivity a in Eq. ([8]) decreases with 
increasing scattering time igc- This can be interpreted as 
a system for which the particles are trapped due to quan- 
tum interference phenomena and conductivity is mainly 
caused by non-elastic collisions, where the collision rate 
increases with decreasing scattering time tsc-— For /? > | 
the reverse behavior can be observed, i.e., the conductiv- 
ity decreases with increasing collision rate. We already 
pointed out that various experiments for real quasicrys- 
tals revealed an electronic transport in agreement to the 
conductivity in the subdiffusive regime (/3 < \)'^— 

The scaling behavior of the mean square displacement 
d{t) is visualized for the silver-mean chain C^s in Fig. 
[5] for s = 1 and different coupling strengths w. The 
results are averaged over 300 initial positions of the wave 
packet around the center of the chain. We observe a 
scaling behavior according to Eq. ^ over several orders 
of magnitude till the width becomes constant due to finite 
size effects. The corresponding scaling exponents /? for 
the average mean square displacement d{t) are shown in 
Fig. [3] for different metallic-mean chains. Thereby, the 
scaling exponents P steadily increase with the coupling 
parameter w and the parameter h of the inflation rule. 

In Fig. [2] the mean square displacement d{t) is also vi- 
sualized for the two- and three-dimensional silver-mean 
labyrinth tiling, where we averaged in two dimensions 
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FIG. 2. Mean square displacement d{t) of a wave packet averaged over difTerent initial positions for the silver-mean chain C^-^ 
(left), the two-dimensional labyrinth tiling (center), and the three-dimensional labyrinth tiling L^'^'^^ (right) for s = 1. 
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FIG. 3. Scaling exponents /3 for the golden-mean (top) and 
the silver-mean (bottom) systems in one, two, and three di- 
mensions for s = 1. All results are averaged over different 
initial positions of the wave packet. 



over 50 and in three dimensions over 20 initial positions 
of the wave packet due to the limitations of computing re- 
sources. The results arc qualitatively very similar to that 
of the quasipcriodic chains and the scaling exponents in 
Fig. [3] even show that the one-dimensional scaling expo- 
nents /3^'^ are nearly identical to the results for the two- 
and three-dimensional labyrinth tiling. 



IV. RG APPROACHES FOR GOLDEN- 
SILVER-MEAN CHAINS 



AND 



The structure and the dynamical properties of the 
metallic-mean systems can be described by an RG ap- 
proach, which is a mathematical apparatus for the sys- 
tematical study of the changes of a physical system 
viewed at different length scales. The general idea is that 
a given system can be transformed (renormalized) in an 
RG step in such a way that we obtain a new system with 
a reduced number of degrees of freedom (i.e. sites and 
bonds). We make use of the fact that the metallic-mean 
chains contain for w — only isolated sites (atoms) and 
biatomic clusters (molecules) coupled by a strong bond as 
visualized in Fig. [4l This yields three highly degenerate 
energy levels: E = ior the atomic sites and E = ±s for 
the bonding and antibonding states of the moleculesi^ 
For non-zero parameters w one finds a coupling between 
these isolated clusters, where the dominant contribution 
occurs between sites of the same type.— For w ^ s also 
the probability density of the wave functions is concen- 
trated either on atomic or on molecular sites depending 
on the energy E. For the Fibonacci chain C^" Niu and 
Nori^^ distinguished two possible RG approaches: 

Atomic RG — The renormalized grid after one RG 



5 



C^"D C«^5 O C^"D C^^J O C^^5 

O C^^^^^^— 

O O 

FIG. 4. Atomic RG (top) and molecular RG (bottom) for 
the Fibonacci chain C*": in the atomic RG all atomic sites 
survive an RG step, and in the molecular RG all molecular 
sites are replaced by a new site. The thick (thin) lines denote 
the strong (weak) bonds. 



step consists of the atomic sites of the original chain as 
shown in Fig. |4l These new sites are then connected 
by new strong and weak bonds. In particular, atoms 
separated by one molecule in the original chain get con- 
nected by a new strong bond and atoms separated by 
two molecules in the original chain become connected by 
a new weak bond (cf. also Fig. [TT]) . 

Molecular RG — In one RG step all molecular sites of 
the original chain are replaced by new atomic sites (cf. 
Fig. SI). Again the atoms in the renormalized grid can be 
connected by new bonds, where the sites of neighboring 
molecules are connected by a new strong bond and sites of 
molecules separated by an atomic site become connected 
by a new weak bond (cf. also Fig. [T2|) . 

Both RG schemes yield a new Fibonacci chain, which 
is scaled in length and energy. Basically, the substitution 
of a cluster of bonds by a new bond during an RG step 
corresponds to the block diagonalization of the Hamilto- 
nian of the original chain. This yields an effective Hamil- 
tonian for the new chain containing the bond strengths 
between the remaining sitesi^ The derivation of the scal- 
ing factors z for the new bond strengths as well as the 
length scalings c under the application of one RG step 
is described in detail in the Appendix lA II for the atomic 
RG and Appendix lA 21 for the molecular RG approach. 

For the silver-mean chain C^s 

we can apply the same 
approach. In the atomic RG the renormalized grid again 
consists of all atomic sites of the original chain (cf. Fig. 
[5]). By appropriately assigning weak and strong bonds we 
obtain a new silver-mean chain, which is scaled in length 
and energy. However, the molecular RG approach for 
the silver-mean chain is more complicated and requires 

CK)CKKK>0<KK><KXKK>00<KK>00<)0-^ 

O C»«D O O C^^D O 

O O 

0<XKK><»<>OOH><X><K>OO<>OCKKK^ 

»D Cm™:^^3 C^^3 O—MD C^^^WD 

O CM~0 O O C»»D O 

FIG. 5. Atomic RG (top) and molecular RG (bottom) for the 
silver- mean chain C^^: in the atomic RG all atoms survive 
and for the molecular RG all molecular sites are replaced by 
a new site, which results in a new silver-mean chain only after 
two steps. 



a two-step process to yield a new silver-mean chain as 
shown in Fig.O The new silver-mean chain contains then 
all atomic sites which are located between two molecular 
sites of the original chain. Details are presented in Ap- 
pendix |X3] for the atomic RG and Appendix lA 41 for the 
molecular RG approach. 



V. RG APPROACH FOR WAVE-PACKET 
DYNAMICS IN METALLIC-MEAN CHAINS 

The scaling exponent /3 of the mean square displace- 
ment d{t) for the golden-mean and the silver-mean chain 
in the regime of strong quasiperiodic modulations {w <C 
s) can be determined by an approach proposed by Abe 
and Hiramoto2ii2^ and further refined by Piechonj^ It is 
based on the idea that the time evolution of an initially 
localized wave packet can be considered as the realiza- 
tions of successive steps of the above described RG. In 
particular, this leads to the following two transforma- 
tions: 

Length scaling — In one RG step the width d(t) of the 
wave packet scales with the grid spacing c, i.e., d{t) — > 
d'it)/c. 

Time scaling — The time evolution of a wave packet is 
governed by the hierarchic structure of the tiling, which 
leads to the occurrence of hierarchic resonances on dif- 
ferent time scales (cf. Fig. [7]).— >2ii^ Since for w <C s the 
energy scales of these resonances are very different, no 
interference of electrons on different time scales is consid- 
ered. Further, the scaling factor for the time t under one 
RG step equals the inverse of the scaling factor z = 1 /t of 
the energy. The second factor is obtained by computing 
the strengths of the new bonds of the different clusters 
in the RG approach by applying Brillouin-Wigner (BW) 
perturbation theory, as discussed in the appendix. 

Therefor, from the scaling behavior d{t) ~ dot^ of the 
mean square displacement according to Eq. one can 
derive the relation^ii^ 

z?-;^. (9) 

in z 

Abe and Hiramoto as well as Piechon presented only 
results for the Fibonacci chain for the two cases of the 
atomic RG and the molecular RG<2i"— In this section 
we also derive an approximation of the scaling expo- 
nent /? for the silver-mean chain and show for the as- 
sociated labyrinth tilings that the corresponding scaling 
exponents approach the one-dimensional results in the 
regime of strong quasiperiodic fluctuations. The analyti- 
cal expressions are also compared with numerical results, 
which show a very good agreement. 

A. Fibonacci chain C^" 

Using the results by Niu and Nori^i^^ (cf. Appendix 
IXT|) . the energy scaling for the weak and the strong clus- 
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ters in the atomic RG is zau = vP' I . With the length 
scahng c5^*°™ = tJ^ they obtained for /3 the expressions^ 



1 t 



oatom ^ 3 In TAu 



(10) 



For the molecular RG the energy scaling for the weak 
and the strong clusters is zau = wjls and the length 
scaling c™°' ~ r^^ . This yields the scaling exponen1j2S 



(Omol 



2 In TAu 
In-^ +ln2 



(11) 



While the results by Abe and Hiramoto are based on a 
quantitative argument, rigorous mathematical consider- 
ations by Damanik22ii2S showed that for a periodic chain 
with a quasiperiodic potential (diagonal model) modu- 
lated according to the Fibonacci sequence a very similar 
relation can be found for u> ^ s. With the constants c\ 
and C2 this result is given by ci / In < /3 < C2 / In . 

In Fig. [6] both scaling exponents, /3au™ a-nd /3au', are 
compared with the numerical results of the scaling expo- 
nent /? of the mean square displacement (i(i). While for 
small coupling parameters w the results are very close 
for both RG approaches, the molecular RG theory yields 
significantly better results for the scaling exponent /3 for 
larger values of w. This can be understood by the fact 
that the dynamics of an initially localized wave packet 
are actually governed by the atomic and the molecular 
RG approachj22iS^ The reason is that independently of 
the original type (atomic or molecular site) of the ini- 
tial position Iq of the wave packet after one RG step the 
new initial site Vq can be either an atomic or a molecular 
site. Thus, the scaling exponent /3 of the mean square 
displacement d(t) shown in Fig. [S] strongly depends on 
the percentages Patom and Pmoi of atoms and molecules 
in the golden-mean chain. As the percentage Pmoi is sig- 
nificantly higher, it is not surprising that the exponent 
/Jau' shows a better agreement with the numerical data 
for the exponent /3. One can also determine a weighted 
average of both scaling exponents, that is, 

/3au — Patom/^Au Pmol /3au (12a) 



AT^^oo PAu j (PAu -PAu ) 



Nr.. 



© TAu - 1 
TAu + 1 



/jatom I 
PAu + 



TAu + 1 



PAu 



(12b) 
(12c) 



A comparison of this expression with the numerical re- 
sults in Fig. [6] shows a very good agreement for w s, 
but for larger w the agreement is not as good as for /^au'- 
In the second case we have to keep in mind that we ap- 
ply perturbation theory and the results are only valid for 
w ^ s. As the RG approach overestimates the results for 
larger values of w, usually the lower of the two exponents, 
/3f ""^ and 13'^°\ yields the better results for w > 0.3s. 

It is also possible to determine the scaling exponents, 
f^Au™ and /3au': numerically. This is achieved by study- 
ing wave packets which are initially localized at a position 
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FIG. 6. Scaling exponent /? of the mean square displacement 
d{t) averaged over different initial positions in comparison 
to the analytical results of the atomic and molecular RG ap- 
proach for the golden-mean (top) and the silver-mean systems 
(bottom) in one, two, and three dimensions for s = 1. 



Iq that survives all atomic/molecular RG steps for an ap- 
proximant of finite size. Using this method Hiramoto and 
Abe found a good correspondence between the theoretical 
and numerical results of the atomic RG for w < s/i^ 



B. Silver- mean chain C^^ 

The derivations for the atomic RG approach in Ap- 
pendix IA3I lead for both types of clusters to an energy 
scaling of ZAg = w/s. With the corresponding length 



scaling Ca^ 



T^g we obtain the scaling exponent 
InTAg 



/Datom 
PAg 



In- 



(13) 



For the molecular RG the energy scaling for the weak 
and strong cluster is different. In Appendix I A 41 we ob- 
tain for the strong cluster = up' j and for the weak 
cluster < j for w <^ s. On average we ex- 
pect to obtain a superposition of the spreading caused 
by both clusters. Using Eq. which relates the num- 
ber of strong and weak bonds/clusters in the silver- mean 
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chain, we can derive an average energy scahng of 

2Ag = ^ (#.(Ca)4'g + #u.{Ca)zTg) 



1 



''Ag wc 



1 + TAg 1 + TAg 



(14a) 
(14b) 



With the length scaling of c™°' 
scaling exponent /3™g' the expression 



Tf,„ we obtain for the 



omol 



2 In TAg 



In (1 + TAg 



In I 



■ '^Ag2;Ag 



(15) 



In Fig. [5] the analytical expressions for the atomic 
and the molecular RG are compared with the numer- 
ical results. For the molecular RG we always plot a 
range of possible scaling exponents P™^ by choosing 
•^Ag S [OjUi^/s^]. The results show that the scaling ex- 
ponent slightly underestimates the numerical values 



for w <C s. In contrast, the analytical values /^^g™ of the 
atomic RG are significantly higher than the numerical 
results. As already discussed for the golden-mean chain, 
the exponent /? depends on the initial position /q of the 
wave packet. For instance. Fig. [7| shows that an atomic 
initial position which remains an atomic site under the 
RG transformation is well described by the atomic RG. 
Hence, on average we expect to obtain a scaling expo- 
nent j3 between the results predicted by the atomic and 
the molecular RG. Using the percentages of atoms and 
molecules in the silver-mean chain according to Eq. ([2|), 
we obtain for the average scaling exponent 
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in correspondence with Eq. (jl2p . This expression is vi- 
sualized in Fig. [6] as well. The results show that the ana- 
lytical values correspond to the numerical values within 
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7. Mean-square displacement d{lo 

->Ag 



4040, t) for the 

silver-mean chain C'^^ for w — Q.l and s = 1. The initial site 
lo of the wave packet remains an atomic site for the first five 
atomic RG steps. The slope is well described by 



the error bars for w <C s. But again, the lower expo- 
nent /^Ag' shows a better agreement for larger w. Fur- 
ther, by comparing the results for the golden-mean and 
the silver-mean chain we find that the exponent ^Ag is 
slightly larger than /3au hi agreement with the numerical 
results in Fig. [31 

Further, the repeating patterns of d{t) in Fig. [7] reflect 
the hierarchic resonances of the wave packet dynamics on 
different time scales. We can also determine the scaling 
factor z directly from the distance of successive patterns. 
In the case of the atomic RG for the silver-mean chain 
for w = 0.1 and s = 1 the patterns are repeated in Fig. 
[7] after a factor of 10 on the time scale, which yields the 
scaling factor ZAg = 1/10 for the energy and, thus, for 
the bond strengths. This result agrees perfectly with the 
theoretical result z = w/s. 



VI. RG THEORY FOR HIGHER DIMENSIONS 

Like the one-dimensional systems, the labyrinth tilings 
possess a hierarchical structure. In Fig. [8] we show the 
corresponding atomic and the molecular RG approach 
for the golden- mean labyrinth tiling In this tiling 

we can identify three different clusters (strong, medium, 
and weak) in correspondence to the three different bond 
types. In general, the clusters of the RG in d dimen- 
sions are given by the product of the corresponding one- 
dimensional clusters due to the product structure of the 
labyrinth tiling (cf. Figs.[T5lto [20| . Hence, the number of 
different clusters depends on the number of bond types, 
which in d dimensions is given by d -I- 1 . 

The derivations of an analytical expression for the scal- 
ing of the length and bonds strengths for the labyrinth 
tiling follow the same way as before. As an example they 
are described for the golden-mean labyrinth tiling in 
the Appendices IB II and IB 21 In two dimensions we obtain 
a grid scaling of c^'°™'^'^ = t^^ for the atomic RG and 



mol,2d 



^Au molecular RG. With the BW per- 

turbation theory we determine the new bond strengths, 
where we obtain a qualitatively different result for the 
atomic and the molecular RG. 

For the atomic RG the scaling factor for the bond 
strength is zj^^ = (zau)^ = w'^/s'^. Plugging this into 
Eq. ([9]) we obtain the same analytical expression as in 
one dimension, i.e., Eq. pO|) . This result is not surprising 
because the scaling factors for the grid spacing and the 
bond strengths arc the squares of the one-dimensional 
results, which cancel each other according to Eq. 
Since the dominant coupling between the edge sites of a 
d-dimensional cluster originates from the bonds along its 
diagonal and the couplings along this diagonal are the 
products of the coupling strengths of d one-dimensional 
clusters (cf. Figs. [15] to [17]), for the atomic RG in d di- 
mensions the scaling factor for the bond strengths and, 
hence, the energies is z'^. 

Although the structure of the clusters is very similar 
for the molecular RG, the result is somewhat different 
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FIG. 8. Strong (black lines), medium (dark gray lines), and weak (light gray lines) clusters for the atomic RG (left) and the 
molecular RG (center) for the golden- mean labyrinth tiling C^^. Sites of the renormalized grid are shown in black and sites 
belonging to more than one cluster in white. In the right panel we show a part of the golden-mean labyrinth tiling which 
after one atomic RG step results in the strong cluster in the left panel. 



than the one-dimensional exponent in Eq. PT|) . The 
derivations in Appendix IB 21 yield a scaling of = 
2 (zAu)^ = u'^/2s^ with an additional factor of 2 originat- 
ing from the normalization condition. Hence, we obtain 
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For very small coupling parameters w — > the first term 
of the denominator dominates and this expression ap- 
proaches the one-dimensional result in Eq. (fTT|) . We plot 
this function in Fig. [H] and find that it is greater than 
the one-dimensional results of the molecular RG for all 
considered values of w. For the average scaling behavior 
we have to consider the contributions of both RG ap- 
proaches according to Eqs. (jTU]) and p7|) . By calculating 
the percentages of atoms and molecules, we obtain the 
average scaling exponent 
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o / oniol,2d /3atom,2d 
I^PAu ^ PAu 



© oatom,2d 
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The result for the three-dimensional labyrinth tiling 
£3d,Au follows in the same way. We already pointed 
out that the atomic RG approach yields the same ex- 
ponent in every dimension (cf. Eq. pO|) ). Due to the 
normalization condition we obtain in the molecular RG 
an additional factor of 4, which yields the scaling fac- 



3d 
Au 



4(2;au) = w^/2s^. Hence, with c 



mol.3d 
Au 



tor Z 

the analytical expression for the molecular RG is 
^moi,3d _ (2lnTAu) / (Hs/w) + ^ln2). This yields the 
average scaling exponent 
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atom, 3d 
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(19) 



In Fig. [5] we compare the numerical and the analyti- 
cal results for the average scaling exponent /3au for the 
golden-mean labyrinth tilings. Within the error bounds 
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FIG. 9. Scaling exponent jS of the mean square displacement 
d{t) averaged over different initial positions in comparison to 
the analytical results of the atomic and molecular RG ap- 
proach for the golden-mean labyrinth tiling (top) and 
the silver- mean labyrinth tiling C^^ (bottom) for s = 1. 
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we find a good correspondence in the regime of strong 
quasiperiodic modulation. Further, although the analyt- 
ical expressions for the scaling exponent /3au in different 
dimensions are quite different, we see that they approach 
each other for w ^ s. This is also in good agreement with 
the numerical data in Fig. [31 which show no significant 
differences in one, two, and three dimensions. 

Also for the silver-mean labyrinth tiling we can 
apply the atomic and molecular RG. Like for the golden- 
mean labyrinth tiling we find three types of clusters in 
two dimensions, which are visualized in Fig. 1101 Again, 
we have to distinguish between the atomic and the molec- 
ular RG approach. For the atomic RG the scaling fac- 
tors of the length as well as the bond strengths are the 
squares of the one-dimensional results, which yields the 
same analytical expression for /3 (cf. Eq. (|13p ) as in one 
dimension. 

For the molecular RG we have to carefully adapt the 
RG approach because the scaling factors of the energies 
are not identical for the different types of clusters. Hence, 
we explicitly have to determine the new average scaling 
factor ZAg (cf. Eq. ^^). The scaling factors for the clus- 
ters are given as the product of the one-dimensional re- 
sults, which yields for the strong, medium, and weak clus- 
ters of the two-dimensional labyrinth tiling the factors 

4g'' = ^v^^ = ^ry/^^ and zif" = {zt/, 

respectively. By considering the percentage of each of 
these clusters in the labyrinth tiling, we obtain the aver- 
age scaling factor 

1 2 

ZAg = {#s{Ca)f ^Ig^'^ + ■T72#s{Ca)#w{Ca)z'^g^'^ 

+ ^i#UCa)fzlf'' (20a) 

According to Eq. © this yields the same scaling expo- 
nent as for the one-dimensional system (cf. Eq. (llSp ). 
The average scaling exponent ^Ag follows then from Eq. 
(jlSp by replacing the quantities for the golden-mean 




FIG. 10. Strong (black lines), medium (dark gray lines), and 
weak (light gray lines) clusters for the atomic RG (left) and 
the molecular RG (right) for the silver-mean labyrinth tiling 
£^^. Sites of the renormalized grid are shown in black. 



labyrinth tiling with that of the silver-mean labyrinth 
tiling. Since we only know the upper bound of the 
scaling factor of the weak cluster in one dimension (i.e. 
-^Ag € [0, if'^/s'^]), the analytical expressions again de- 
scribe a range of possible exponents. A comparison with 
the numerical results in Fig. [H] shows that the final scal- 
ing exponent is very close to that of the atomic RG 
and, hence, overestimates the actual results. Rather, the 
molecular RG works better than the average. By per- 
forming the same calculations also for three dimensions, 
we obtain that the scaling exponent /3^g is even closer to 
the scaling exponent of the atomic RG. Hence, we would 
expect that also the numerical results for /? approach the 
analytical expression for /3^'g ™ in Eq. ([T3|) with increas- 
ing dimensionality. However, we do not observe such 
a behavior. A possible explanation for these discrepan- 
cies could be that determining an average scaling factor 
for the bond strength in the molecular RG according to 
Eq. (PO]) from the energy scalings of the different clus- 
ter types does not adequately display the dynamical be- 
havior. At least, the numerical and analytical results 
approach each other with increasing quasiperiodic mod- 
ulation, and we find a correspondence for w = 0.1 within 
the error bounds. 



VII. CONCLUSION 

We have studied the quantum diffusion in quasiperi- 
odic systems in one, two, and three dimensions by in- 
vestigating the time evolution of wave packets. By nu- 
merical calculations of the scaling behavior of the mean 
square displacement d(t) of a wave packet we have ob- 
served the occurrence of anomalous transport for all cou- 
pling strengths w S (0, s) and that the scaling exponents 
/? are more or less independent of the dimension. 

Further, we have extended an RG method, originally 
proposed by Abe and Hiramoto for the golden-mean 
chain to obtain an analytical expression also for the 
scaling exponent /3 of the silver-mean chain. With the 
same method we have been able to show that the scal- 
ing exponents f3'^'^ of the labyrinth tilings approach the 
one-dimensional scaling exponents /S^'^ for the golden- 
mean and the silver-mean system in the regime of strong 
quasiperiodic modulation {w <C s). The analytical re- 
sults are also in good agreement to our numerical results. 

Although properties like the structure of the energy 
spectrum strongly depend on the dimensionality, these 
derivations also showed that the product structure is an 
essential reason for the rather similar scaling exponents 
(3 in different dimensions. Hence, further research of the 
transport properties is needed especially with a focus on 
non-separable tilings. Nevertheless, we have been able to 
draw a connection between the structure of the quasiperi- 
odic systems and their transport properties with the RG 
approach. We found that for strong quasiperiodic mod- 
ulations the wave-packet dynamics are governed by the 
hierarchical structure of the quasiperiodic systems lead- 
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ing to the occurrence of hierarchic resonances on different 
time scales. 

AdditionaUy, for w <C s the numerical and analytical 
results for the exponent /3 show the occurrence of sub- 
diffusive wave-packet dynamics (/? < 1/2). According to 
the generalized Drude equation this leads to a decrease 
of the conductivity with increasing structural order as 
it is observed in real quasicrystalsii^ Hence, labyrinth 
tilings can be useful models to study the characteristics 
of higher-dimensional quasiperiodic systems efficiently by 
numerical methods and also analytically. 



Appendix A: Renormalization of lengths and 
energies of quasiperiodic systems 

This section comprises the results for the sealing of 
lengths and energies in the RG approach introduced by 
Niu and Norii^i^^ Although the scaling factors have been 
derived for the Fibonacci chain before ;^^'^^ we will briefly 
sketch the calculation for the atomic RG and mention 
the results for the molecular RG because we make use of 
them also for the octonacco chain. Further, it allows us 
to compare the results for different dimensions. 



Fibonacci chain C 



Atomic RG 



In the atomic RG approach for the Fibonacci chain 
the renormalized grid after one RG step consists of the 
atomic sites of the original chain (cf. Fig. [4]), which leads 
to a scaling of the bonds according to the substitution 
rule in Fig. [TTJ The scaling of the grid spacing is given 
by the substitution matrij*^ of the corresponding RG ex- 
pansion ("Paai"™)"^' which relates the number of symbols 
w and s in successive RG steps: 
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FIG. 11. Substitution rule Pau" of the atomic RG for the Fi- 
bonacci chain C*" shown in Fig. |4] (top) with the energy scal- 
ing for the strong cluster (left) and the weak cluster (right). 



The Perron-Frobenius eigenvalue Apr of this matrix cor- 
responds to the inverse of the scaling factor c^*°™ of the 
grid spacing, i.e.. 



det(R - AI) 



atom 



A2- 

1 
Apr 



4A - 1 = 
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V5 
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The scaling of the new bond strengths is calculated by 
Brillouin-Wigner (BW) perturbation theory for each of 
the two substitution clusters in Fig. [11] For the metallic- 
mean systems the coupling parameter w is treated as a 
perturbation. Hence, the Hamiltonian H = Hq -I- Hi is 
decomposed into a Hamiltonian Ho(w — 0) of the unper- 
turbed system and a perturbation Hi = H — flo{w = 0). 
The approach is described in detail by Niu and Nori)^ 



a. Energy scaling for strong cluster 



The Hamiltonian of the strong cluster in Fig. [TTjis 
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Solving the Schrodinger equation for the unperturbed 
Hamiltonian Hp yields the wave functions |^'i) = |1), 
14-2) = |4), and 1^-^) = ^ (|2) ± |3)) with their corre- 
sponding energy values i?i/2 = 0, and = ±s. 

The strength s' of the new bond is given by the ma- 
trix element (5'i| H |^'2), where only the leading non-zero 
term of the perturbation expansion is relevant. The BW 
perturbation theory yields no contribution in the first 
order because (^'i| Hi 1^*2) = 0. In the second order we 
obtain a new coupling strength of 
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Here, P denotes the projection operator out of the sub- 
space for a given g-times degenerate eigenstate 5'^ with 
the energy of Hq, i.e., for the strong cluster it is given 



byP 



in 



6. Energy scaling for weak cluster 

Analogously this approach is applied to the weak clus- 
ter. The Hamiltonian Hq yields the atomic wave func- 
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tions l^-i) = |1), 
as well as the 
^(|2)±|3)) and 

±s. 



'^2) = |6) with the energy E1/2 
molecular wave functions |^^) 



.4=) = ^ (|4) ± |5)) with the ener- 
gies Ef/^ 

The strength w' of the new bond is again given by 
the matrix element H 1^2) • We obtain no contribu- 
tion in the first order expansion due to (^"11 Hi [^'2) ~ 



and with the projection operator P 
der term vanishes, HiP 



m){n\ + 

\'i>A ) (4'/ I also the second or- 



E-U, 



■PHi l^-a) = 0. In the 



third order we have four non-zero terms: 
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Hence, for the atomic RG of the Fibonacci chain the 
scaling of the bond strength of both clusters and, thus, 
the scaling of the energy is given by zau = j . 



Fibonacci chain C 



Molecular RG 



In the molecular RG of the Fibonacci chain all molecu- 
lar sites of the original chain are replaced by new atomic 
sites as shown in Fig. IH which leads to a scaling of the 
bonds according to the substitution rule 7^™°' i'^ I^iS- GSI 
The scaling factor c^u' of the grid spacing is again deter- 
mined by the substitution matrix 
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FIG. 12. Substitution rule Vtt of the molecular RG for the 
Fibonacci chain C^" in Fig. |4] (top) with the energy scaling 
for the strong cluster (left) and the weak cluster (right). 



3. Octonacci chain — Atomic RG 

Like in the atomic RG approach for the Fibonacci 
chain, here only atomic sites survive during an RG step 
as shown in Fig. [5l The corresponding substitution rule 
and the scaling of the clusters are shown in Fig. [131 The 
substitution matrix is now given by 



atom\ l 
Ag ) 
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and thus the grid spacing scales with 
1 1 
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The bond strength for the strong cluster is very simple, 
and for the weak cluster (cf. Fig. [T3|) we already derived 
its strength in Eq. (|A4[) while studying the strong cluster 
in the atomic RG of the Fibonacci chain (cf. Fig. [TT|) . 
i.e., 
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The scaling factors for the new bonds of the two sub- 
stitution clusters in Figs. [12] follow by BW perturbation 
theory. The derivations are described in detail by Niu 
und Nori^^ and yield 
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This corresponds to an energy scaling zau — w /2s. 
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FIG. 13. Substitution rule Vf""^ of the atomic RG for the 
silver-mean chain C^^ in Fig. [5] (top) with the energy scaling 
for the strong cluster (left) and the weak cluster (right). 
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Hence, we obtain an energy scaling ZAg = w/ s. 



4. Octonacci chain C^^ — Molecular RG 



Applying the molecular RG to the silver-mean chain 
C^g leads to some challenges because two RG steps are 
needed to obtain a new silver-mean chain (cf. Fig. [5]). The 
substitution matrix of the corresponding RG expansion 
(cf. Fig. [14]) is given by 
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2 1 



and we obtain a scaling of the grid spacing of 
1 1 



mol _ 



Apr 3 + 2^/2 
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The scaling of the bond strengths is again calculated us- 
ing the BW perturbation theory for each of the two sub- 
stitution clusters in Fig. [HI 



a. Energy scaling for strong cluster 

Here, we have the same substitution rule as for the 
strong cluster in the atomic RG theory for the Fibonacci 
chain (cf. Fig. [TTjl . that is, we get 



w w 

s 



(AM) 



This results in an energy scaling of = j for the 
strong cluster. 



6. Energy scaling for weak cluster 



sites of the cluster and some of the inner sites belong 
to the same subspace. Hence, we have to resolve this 
degeneracy first to compute the new Hamiltonian matrix 
elements. However, the new set of wave functions shows 
already a coupling of the edge sites of the cluster so that 
we do not obtain information about the energy scaling of 
the bonds. 

However, we can obtain an upper bound replacing 
the central weak bond of the weak cluster by a strong 
bond. For this modified cluster we can apply the RG 
approach. In particular, we obtain for the Hamilto- 
nian Hp of the modified cluster the atomic eigenstates 
= |1), = |8) with Ei/2 = and the molecular 
eigenstates |vI/±) = ^ (|2) ± |3)), = (|4) ± |5)), 

|*±) = ^(|6)±|7))withi?3±_5 = ±s. 

The leading non-zero contribution occurs only in the 
fourth order of the perturbation expansion: 
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This new cluster possesses a larger coupling between 
the edge states than the original cluster shown in Fig. 
1141 Hence, the corresponding scaling factor provides an 
upper bound for the scaling factor of the weak cluster 
with < w'^ /s^. 



The computation of the strength of the new weak bond 
turns out to be rather complicated. The reason is that 
we cannot apply the method used so far because the edge 



Appendix B: Renormalization of Lengths and 
Energies for the Golden-Mean Labyrinth 
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FIG. 14. Substitution rule P^g of the molecular RG for the 
silver-mean chain C^^ in Fig. [5] (top) with the energy scaling 
for the strong cluster (left) and the weak cluster (right). 



1. Atomic RG 

For the labyrinth tiling we use the same technique as 
for the quasiperiodic chains. For the atomic RG the 
renormalized labyrinth tiling contains only the atomic 
sites of the original tiling as shown in Fig. [8] for the 
golden- mean labyrinth tiling C^^. In two dimensions we 
obtain three types of clusters. The corresponding sub- 
stitution rules for the strong, medium, and weak cluster 
are shown in Figs. [151 [HI and [iTl The scaling of the 
grid spacing c can be easily obtained from the scaling 
N'g^ — 7> cNa of the one-dimensional system size. Thus, we 
obtain in d dimensions 
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Hence, for the atomic RG in two dimensions the scahng 
of the grid spacing is 



atom. 2d 



(B2) 



a. Energy scaling for strong cluster 



We assign each vertex {l,m) of the two-dimensional 
clusters a unique index r = [|] + [(m— 1)^], where / 
is the length of the underlying one-dimensional cluster. 
This corresponds to the row-wise assignment of succes- 
sive indices beginning from the lower left corner to the 
upper right corner. With this mapping we can write the 
wave functions as a vector and the Hamiltonian of the 
strong cluster of the RG is given by 
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The corresponding Hamiltonian Hq yields the normalized 
wave functions = |1), |*2) ^ |2), 1*3) = |4), 1*4) = 
|5), \^,) = |7), 1*6) = |8), and 1*^) = ^ (|3) ± |6)) 

with the energy values i?i-6 = and Ef — ±s^. 

The strength of the new bond is given by the matrix 
element (*i| H l^'e)- In analogy to the one-dimensional 
case, the first order contribution is (vPil Hi j^I/g) = 0, 
and the second order contribution yields with P = 
1*7 ) (*7 I + 1*7 > (*f I the non-zero term 
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FIG. 15. Substitution rule for the strong cluster in the atomic 
RG approach of the golden-mean labyrinth tiling 
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FIG. 16. Substitution rule for the medium cluster in the 
atomic RG approach of the golden-mean labyrinth tiling £^". 



b. Energy scaling for medium cluster 

For the cluster of medium strength we use the same 
approach. The Hamiltonian again contains all couplings 
as shown in Fig. [TBI and the eigenstates of Hq are |*i) = 
|1), 1*2) = |2), 1*3) - |3), 1*4) - |6), 1*5) = |7), 1*6) = 
|10), 1*7) = 111), 1*8) = |12) with Ei_s = and |*^j 
^ (|4) ± |8)), |*± ) = ^ (|5) ± |9)) with E^^^^ . 

The new bond strength of the medium cluster is given 
by the matrix clement (*i|H|*8). Thus, there is no 
contribution in the first and the second order expansion. 
In the third order we have four contributing terms: 
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c. Energy scaling for weak cluster 

For the weak cluster the Hamiltonian again contains 
all coupling strengths as shown in Fig. [ITl The Hamilto- 
nian Ho yields the atomic wave functions |*i) = |1), 
1*2) - |2), 1*3) ^ |3), 1*4) = |6), 1*5) = |7), 
1*6) = |12), 1*7) = |13), 1*8) = |16), 1*9) = |17), and 
1*10) = 1 18) with an energy value i?i_io = as well 
as the molecular wave functions |*n) = (|4) ± |8)), 

|*±) = -^{\5)±\9}), 1*^3) = :i5(|10)±|14)), and 



V2 



V2 



l*w) = 71 (111) ± |15)) with energies E^,_,^ = ±s\ 

The strength of the new bond for the weak cluster is 
given by the matrix element (*i|H|*io). As for the 
medium cluster we obtain no contribution in the first 
and the second order expansion. In the third order we 
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FIG. 17. Substitution rule for the weak cluster in the atomic 
RG approach of the golden-mean labyrinth tiling 



have four contributing terms: 
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Hence, for the atomic RG of the golden-mean labyrinth 
tiling the scaling of the bond strengths for all three clus- 



ters is given by zj^^ = (zau) 



2. Molecular RG 

Also in the molecular RG of the golden- mean labyrinth 
tiling three types of clusters occur (cf. Fig. [8]). The corre- 
sponding substitution rules for the strong, medium, and 
weak clusters are shown in Figs. [THl [121 and [201 respec- 
tively. Again, the length scaling in two dimensions fol- 
lows from Eq. (|Bip as 
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a. Energy scaling for strong cluster 
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FIG. 18. Substitution rule for the strong cluster in the molec- 
ular RG approach for the golden-mean labyrinth tiling 



{'i'll Hi l^'J) = 12. Hence, the scaling of the energies 



l\2 



of this cluster is given by (s') 



2^ ^2 



h. Energy scaling for medium cluster 

The Hamiltonian for the medium cluster contains all 
coupling strengths as shown in Fig. [191 The eigenstates 
of Ho are l^i) = |2), \^2) = |7) with E1/2 = and 



^(|1>±|4)), 



^(13) ±15)), 



^ (|6) ± |9)), |*±) = ^ (|8) ± |10)) with Et_, = ±s^. 

The strength of the new bond for the bonding state 
IS given by the matrix element (*;^|H|*^). The BW 
perturbation theory yields no first order contribution be- 
cause I Hi I'^f) = and the second order contribu- 
tion for the bonding state is given by 
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c. Energy scaling for weak cluster 

For the weak cluster the Hamiltonian again contains 
all coupling strengths as shown in Fig. [20l The Hamilto- 
nian Ho yields the atomic states |5'i) = |2), 1^*2) = |6), 



The Hamiltonian of the strong cluster of the molecular 
RG contains all coupling strengths shown in Fig. [181 and 
the Hamiltonian Hq yields the four molecular eigenstates 
\^f) = ^(|1)±|3)), = ^(|2)±|4)), ^ 

^ (|5) ± |7)), and |*±) = ^ (|6) ± |8)) with i?± = ±s\ 
The new bond strength for the bonding state is 
given by the matrix element (^f|H|\['^), where only 
the leading non-zero term of the perturbation expan- 
sion is relevant. The BW perturbation theory already 
yields a contribution in the first order expansion with 
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FIG. 19. Substitution rule for the medium cluster in the 
molecular RG for the golden- mean labyrinth tiling £^". 
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l^fg) = |7), 1*4) = |8) and {^5) = |12) with E = as well 
as the molecular wave functions l^*^) = (|1) ± |4)), 

. ^(|3)±|5)), = ^(|9)±|11)), and 

-1= (|10) ± |13)) with the energies ±5^. 

The new bond strength for the bonding state follows 
from the leading non-zero term of the perturbation ex- 
pansion of (vEfjl"! H The BW perturbation theory 



yields no first order contribution Hi I'i'g) 

the second order contribution is given by 
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Hence, the molecular RG yields for all three types of 
clusters the scaling factor z"^^ = w^/2s^. 
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FIG. 20. Substitution rule for the weak cluster in the molec- 
ular RG approach for the golden-mean labyrinth tiling 
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